Effect of anisotropic band curvature on carrier multiplication in graphene 
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We study relaxation of an excited electron in the conduction band of intrinsic graphene at zero 
temperature due to production of interband electron-hole pairs by Coulomb interaction. The elec- 
tronic band curvature, being anisotropic because of trigonal warping, is shown to suppress relaxation 
for a range of directions of the initial electron momentum. For other directions, relaxation is allowed 
only if the curvature exceeds a finite critical value; otherwise, a nondecaying quasiparticle state is 
found to exist. 



I. INTRODUCTON 

Carrier multiplication is a process in which a sin- 
gle photon, absorbed by a material, produces several 
electron-hole (e-h) pairs. Typically, this happens when 
the primary photoexcited e-h pair produces a number of 
secondary pairs of smaller energy via electron-electron 
collisions. This process is very important for optoelec- 
tronic applications: the more e-h pairs are produced by 
a single photon, the more efficiently one can convert light 
into electric current. Graphene is an obvious candidate 
for efficient carrier multiplication, since (i) it has wide 
electronic bands and no energy gap, and (ii) electron- 
electron scattering can be much faster than electron- 
phonon scattering. Indeed, the dynamics of photoexcited 
carriers in graphene has become a subject of many stud- 
ies, both theoretical and experimental.'^ES 

In spite of the numerous studies, this dynamics is still 
not fully understood. Notably, the most basic issue, that 
of the role of electron-electron collisions in the relaxation 
of a single photoexcited carrier in the intrinsic graphene, 
is still under debate. It is well known (see, e. g., Ref. 
that due to simultaneous energy and momentum conser- 
vation, decay of quasiparticles is allowed or forbidden, de- 
pending on the curvature of the quasiparticle spectrum. 
In the context of graphene, the Dirac spectrum is linear, 
so it is exactly on the borderline between the two cases.f^Sl 
Indeed, if an electron in the conduction band with mo- 
mentum p and energy v\p\ {v being the Dirac velocity) 
is scattered into the state with another momentum p' 
and lower energy f |p'| < ^IpI, creating a hole in the va- 
lence band with momentum ph and another electron in 
the conduction band with momentum Pe, the momentum 
and energy conservation conditions, 

P = p'+Pe + P/i, (la) 
v\p\^v\p'\+v\p,\+v\ph\, (lb) 

are compatible only in the special case when all vectors 
lie on the same line. Differe nt ways to r esolve the uncer- 
tainty have been advocated.lilifiliallllllMl Note that the 
collision process in question is precisely the one respon- 
sible for carrier multiplication (Fig. [l]) . 

It might seem that importance of the above-mentioned 
problem of a single excited carrier is limited to low pho- 
toexcitation intensities. However, if many e-h pairs are 
created under intense photoexcitation, it is important to 
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FIG. 1. Various collision processes for an excited electron in 
the conduction band. The straight lines represent the Dirac 
cones and the circles represent electrons. Processes (a), (b) 
conserve the number of carriers in the conduction and va- 
lence bands separately. Processes (a)-(c) are forbidden by 
the Pauli principle in the intrinsic graphene at zero temper- 
ature. Process (d) corresponds to carrier multiplication (the 
initial electron produces one more electron and one hole). 



know whether the population equilibration between the 
valence and the conduction band happens as quickly as 
the thermalization within each band. For example, only 
if the interband population relaxation is slow enough, a 
population inversion between the bands can be achieved, 
and one can think about lasing. Among various two- 
electron collision processes, shown in Fig.[TJ only the pro- 
cesses (d) (carrier multiplication) and (c) (Auger recom- 
bination, reciprocal to the multiplication) can transfer 
electrons between the bands. If these are suppressed, 
three-particle collisions are required to equilibrate the 
populations of the conduction and the valence bands. 

In the present work, we study relaxation of an excited 
electron in the conduction band of intrinsic graphene at 
zero temperature due to electron-electron collisions (the 
process (d) in Fig. [T]), going beyond the Dirac approxi- 
mation and taking into account the electronic band cur- 
vature. Because the curvature is anisotropic due to the 
trigonal warping, the result turns out to depend on the 
direction of the initial electron momentum, as was also 
noted in Ref. |29l For a certain range of directions, the 
process is forbidden. For the directions when relaxation 
is allowed, we calculate its rate. 

If the curvature is weak, the problem corresponds to 
that of a discrete state coupled to a continuum whose 
density of states is abruptly cut off precisely at the en- 
ergy of the discrete level, which can be viewed as a special 



2 



case of the Fano problem!^ Indeed, in the sector with the 
fixed total momentum p, the single-particle excitation is 
a discrete level with the energy vp. The three-particle 
density of states vanishes below this energy, and exhibits 
a steplike discontinuity at the energy vp. In this situa- 
tion, the presumable decay of the discrete state into the 
continuum cannot be described by the Fermi golden rule, 
since the latter is valid only when the density of states 
in the continuum is a smooth function of energy at the 
position of the discrete level. Below it is shown that 
the quantum-mechanical level repulsion between the dis- 
crete state and the continuum plays a major role in this 
problem. As a result, a nondecaying quasiparticle state 
exists (that is, its lifetime being determined by mecha- 
nisms other that electron-electron interaction). Still, the 
quasiparticle spectral weight is reduced, some part of it 
being transferred to the continuum of the multiparticle 
excitations. The quasiparticle state can relax by pro- 
ducing electron-hole pairs only when the band curvature 
along the allowed directions exceeds a finite critical value, 
needed to overcome the level repulsion. 



II. CALCULATION 

To derive the results, outlined above, let us describe 
the electrons by a two-component column fermionic field 
operator V'a(r), where a — 1,...,N labels electronic 
species (the valley and spin degeneracy in graphene cor- 
respond to iV = 4). The Hamiltonian is given by 

N 



p{r) /5(r') 



(2) 



e 
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The first term in Eq. ^ represents the kinetic energy of 
electrons, described by the single-particle Hamiltonian, 

CoP^, 
(3) 

(7a;, ay being the Pauli matrices. It determines the single- 
particle dispersion relation to 0{p^), which we write as 

Ep ^VPT CoP^ T CsP^ cos 3ip = vp± C^p^ , (4) 

where ep(ep) is the energy of an electron in the conduc- 
tion band (a hole in the valence band) with momentum p, 
and (fi is the polar angle of p. The coefficients can be esti- 
mated from the tight-binding modeP^^: the Dirac velocity 
V ~ 7eV- A, the trigonal warping coefficient ^3 w 5eV- A , 
and the electron-hole asymmetry Co ~ 0-8 eV • A . Ne- 
glecting the terms proportional to Co i Ca in Eq. (|3| corre- 
sponds to the Dirac approximation h(p) — vp ■ cr. 

The last term in Eq. ([2| describes Coulomb interaction 
between the electrons, with the electronic density 



N 



P(r) = 



(5) 



and the background dielectric constant Sb of the sub- 
strate incorporated into e^. The dimensionless Coulomb 
coupling strength e^/w can be small if e?, is large enough. 
The largest value of e^/u ~ 2.2 is attained for a graphene 
sheet suspended in vacuum, £(, = 1 . 

The quasiparticle decay rate is given by — 2I]p(e), 
where I]p(e), the retarded self-energy projected on the 
eigenstate of the single-particle Hamiltonian with mo- 
mentum p, should be taken at e corresponding to the 
quasiparticle pole of the Green's function. In the first ap- 
proximation, it can be taken on the mass shell, e — or 

e = — fi-p- The lowest order of the perturbation theory in 
the Coulomb interaction, which contributes to ImEp(e), 
is the second one. It describes decay of the one-particle 
state (an electron or a hole) into three-particle excitations 
(an electron or a hole plus an e-h pair). It turns out, how- 
ever, that in the Dirac approximation, the second-order 
Iml]p(e) has a steplike discontinuity on the mass shell, 
invalidating the simple Fermi golden rule recipe for cal- 
culating the decay i&te^ Explicitly, at |e — t;p| ^ vp, 



Imi;p(e) = -7p6'(e - vp), 



7p 



.N 2 
-'24-^"^ 



— vp, 

V 



(6a) 
(6b) 



9{x) being the step function. In Eq. (6b), A^/24 



comes from the bubble diagram (the direct term), while 
ln(2/e^/'^) comes from the exchange diagram (see Ap- 
pendix |A^ for details). At A^ = 4, the exchange term is 
more than six times smaller than the direct term. 

One may consider higher orders of the perturbation 
theory in /v, while remaining within the Dirac approx- 
imation. Decay into {2n + l)-particle excitations, which 
involves three-particle excitations as virtual intermediate 
states, is suppressed by energy and momentum conserva- 
tion as cx {e — vp)"'9{e — vp) (see Appendix A 2 1, so inclu- 
sion of many-particle excitations does not help to resolve 
the uncertainty. Dressing the decay into three-particle 
states by higher-order corrections can be performed by 
treating 1/N as a formal small parameter. This selects 
the random-phase-approximation (RPA) sequence as the 
dominant subcl ass of d iagrams. In RPA, Im Sp( e = v p) 
strictly vanishe^^IEIISni^ Exphcitly (see Appendix A 6 1 , 



IinSp(e) 



4(e — vp) 

TTiV 



9(e — vp) In 



Crpa 
e ~ vp' 



(7) 



where the upper cutoff of the logarithm is ^rpa = 
mhi{vp, {Ne^ /v)'^vp}, and Eq. Q is valid at e — <C 
^RPA. If Ne^/v <^ 1, the perturbative expression of 
Eqs. ( [6a| ), (6b I is still valid in the parametric region of 
energies {N^Jv)'^vp <C e — vp <C vp, while for Ne^ /v ^ I 
Eqs. (6a I, (6b) are never valid, and only Eq. ([t]) holds. 

Beyond the Dirac approximation, we take into account 
the terms proportional to p^ in Eq. ([s]). Assuming them 
to be small compared to the main Dirac term vp, we 
neglect them wherever they produce small corrections to 
regular expressions (e. g., corrections to the eigenstates 



-0.5 





e-vp 



0.5 



-3 



-1 





at-vq 



FIG. 2. Smearing of the steplike discontinuity in the second- 
order self-energy for (^p > at = 4. The dashed and 
solid curves represent — ImEp(e) in and beyond the Dirac 
approximation, respectively. For < the curve should be 
rotated by 180°. 



FIG. 3. Smearing of the square root singularity at lj — )■ in 
the imaginary part of the polarization operator, — Imllq(Lij). 
The dashed curve represents the Dirac approximation which 
gives l/\/Z) th e so lid curve represents the function ImP(z), 
defined in Eq. (|9bl. 



of the single-particle Haniiltonian), and take them into 
account only in those terms vifhich are singular at e — > vp. 
In the second order of the perturbation theory, instead 

of Eq. we obtain (see Appendix lA3] for details of \Cv\P^- We give its explicit value on the mass shell only. 



where E{m) and K{m) are the complete elliptic integrals. 
The function V{z) is plotted in Fig. |3] Finally, Iml]p(e) 
is evaluated in Appendix |A 7| It vanishes for e < vp — 



the calculation) 

Ccp>0 : Im Ep(e) = -7p 



N 



e ~ vp 



C<p < : Iml]p(e) = -7p 



1 - I - 



e — vp 



(8) 



where the function Tjq{z) is defined in Appendix A3 
The result of its numerical evaluation for TV = 4 is plotted 
in Fig. [2] J^n{z) = for z < — 1 . On the mass shell, 
e — vp + C(pP^j -Fn ~ 1- So, electronic relaxation is 
allowed if cos3iy9 < — C3/C0 ~ 0.16, with the rate 27p, 
and forbidden for cos3(y9 > — Ca/Co- For a hole in the 
valence band with momentum p, the conditions are just 
the opposite _ 

) is vahd for {Ne^ /v^vp <t: 
■^vw ~ {Ne^/vYvp and iV > 1, the 
self-energy should be calculated from RPA in the pres- 
ence of the p^ terms. For this, one first has to calculate 
the polarization operator nq(w) (the effect of non-Dirac 
dispersion on the polarization operator in graphene was 
also studied in Refs. [34l and l35|) . Here we take the p^ 
terms into account only near the singularity at w vq, 
which becomes smeared as (see Appendix A 5 for details 
of the calculation) 



Just like Eq. (6a), Eq. 
e — vp'tivp. When | (^p \p'^ 



nq(w Kvq) = - 



1 



Vl^\, (9a) 



V{z) 



37r 



il + z)K 



1 - z 



- zE 



I- z 



(9b) 



which has a compact form: 



ImEp(e;) = - 



64 



(10) 



The above results [Eqs. (Is]) and ( 10 1] correspond to 
the Fermi golden gule (FGR) with lowest-order or RPA- 
dressed transition matrix elements. They are valid in the 
case ICipb^ ^ iniii{7p, wp/A^}, when the singularity in 
Iml]p(e) is strongly smeared by the band curvature. Let 
us study the opposite case, when the dominant energy 
scale is 7p itself, e. g., near the directions where C^j = 0. It 
should be recalled that FGR works only when the density 
of the final states of the decay (three-particle excitations 
in the present case) is approximately constant, in which 
case the quasiparticle spectral peak has the Lorentzian 
shape. When the density of states is not smooth, FGR- 
based approaches^ ^^'^^ are not valid, and the quasi- 
particle spectral peak is manifestly non-Lorentzian. To 
determine the quasiparticle properties in the non-FGR 
regime, let us study the single-particle (retarded) Green's 
function G'p(e), which determines the quasiparticle spec- 
tral function, — (I/tt) ImG'p(e). 

Let us first analyze the most "dramatic" case when 
Iml]p(e) is given by Eq. (6a). In this case, the retarded 



Green's function (or, more precisely, its projection on the 
eigenstate of the single-particle Hamiltonian with mo- 
mentum p) is given by 



Gp(e) 



e — vp 



7p 



In- 



max 

- vp\ 



ilpO{e - vp) 



(11) 
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FIG. 4. The quasiparticle spectral function (in the units 
of 7p), as determined from Eq. (Ill for ^max/7p = 10. 



The vertical peak represents the (5-function with the spectral 
weight Zp ^ 0.72. 



with the real part of the self-energy reconstructed from 
the Kramers-Kronig relation, and ^max ~ vp determines 
the ultra violet cutoff of the logarithmic divergence, since 
Eq. (6a I is valid for e — vp vp onljB^. G'p(e) has a 
real pole at e — vp = — Ap < 0. Its existence immedi- 
ately follows from the fact that ReEp(e) < 0, which, in 
turn, is a consequence of the usual quantum-mechanical 
level repulsion: the quasiparticle level is repelled from 
the three-particle continuum. Note that introduction of 
any infinitesimal broadening of the step function does not 
affect this result at all. The real pole corresponds to a 
quasiparticle state with an infinite lifetime. Even though 
the quasiparticle does not decay into the continuum, the 
latter still takes away part of the spectral weight, which 
manifests itself in the residue Zp < 1 at the pole. With 
logarithmic precision we can evaluate 



Tp 1 ^max 

— In , 

^ 7p 



1 

l + l/ln(e,„ax/7p)' 



(12) 



In RPA, suppression of ImSp(e), Eq. ([t]), at e — up ^ 
Crpa cuts off the logarithmic divergence in ReS]p(e 
vp) at e~vp-^ Crpa [here ^rpa = mm{vp, {Ne^/vYvp} 
is the same as in Eq. ([7|]. Since ^rpa ^ 7p, we can 
take Re Sp (e = vp) to find the quasiparticle pole, Ap « 



Re Sp(e — vp): 



A. 



vp 



Crpa 

(iVeV-f > 1) 



{Ne^v < 1), 



(13a) 
(13b) 



where Cmax/^RPA ~ {Ne^lv)~\ When Ne^/v < 1, the 
logarithm is large, so E q. (|13a| ) has logarithmic precision. 
When Ne^ /v ^ 1, Eq. ( |13b[ ) represents just an order-of- 
magnitude estimate obtained by plugging Eq. ([7| into the 
Kramers-Kronig relation, and integrating from e — vp = Q 
to e — vp ^ vp. To find the residue at the pole, we note 
that at < vp~e ^ Crpa, from Eq. ([t]) and the Kramers- 
Kronig relation one obtains 

5rpa 

9S]p(e) _ 4 f ^H^RPA/Od^ ^ 2 Crpa 



de TT^N J {£, + vp ~ e)^ n'^N vp — 



(14) 
(15) 



which gives 



(We remind that Ec^. (11 1 is valid only when 7p/^max 
iV(e^/u)^ <C 1, so that the logarithm is large). The spec- 
tral function for ^max/7p = 10 is plotted in Fig. |4] 

If now one gradually increases the band curvature 
C^p^ > 0, (i) the bare quasiparticle level is shifted to for < 0. Consequently, the requirement for the band 



with logarithmic precision. The critical value of the 
band curvature, when the quasiparticle level enters the 
continuum and acquires a finite decay rate is given by 
C^P^^Ap/2. 

It is important that the existence of the infinitely nar- 
row quasiparticle peak in the Dirac approximation, ob- 
tained above using the perturbation theory in e^/v and 
1/N expansion, is, in fact, more general that these ap- 
proximations. Indeed, the existence of the peak follows 
from two facts: (i) Im Sp(e < vp) = 0, which holds in any 
order of the perturbation theory because of energy and 
momentum conservation as discussed in Appendix |A 2[ 
and (ii) Re Sp(e — >■ vp) < 0, which is a consequence of the 
level repulsion between the single-particle and the three- 
particle states. When the band curvature is included, 
the quasiparticle peak and the continuum are pushed to- 
wards each other for > 0, and away from each other 



e — vp + CfP^ ' (ii) the spectral boundary of the contin- 
uum is shifted to e = vp — CipP'^i s-nd (iii) the logarith- 
mic divergence in ReEp(e — > vp) at e ^ vp is cut off 
at e — vp C^p'^. The dressed quasiparticle level en- 
ters the continuum at some critical value of C^p^, needed 
to overcome the level repulsion. This critical value can 
be determined from the condition G~^{vp + CvP^) — 0, 
and with logarithmic precision, it is given by C^tpP^ — 
[7p/(27r)] ln(^max/7p)- At this point the quasiparticle 
state acquires a finite decay rate, whose value for suf- 
ficiently large CwP^ is given by — 2ImI]p(ep), as deter- 
mined by Eq. 



curvature to exceed a finite critical value in order to over- 
come the level repulsion and to produce quasiparticle de- 
cay, is also more general than the approximations used 
here. The calculation of critical value itself, of course, 
does rely on approximations. 

Still, one cannot exclude the appearance of nonzero 
Iml]p(e < vp) in the Dirac approximation due to non- 
perturbative effects. For example, nonperturbative gen- 
eration of spectral weight in Imnq(Lij) at uj < vq due 
to excitonic effects has been discussed in Ref. |37l even 
though the validity of these results has been questioned 
in Ref. jiS; This issue calls for further investigation. 




FIG. 5. Directions of the electron momentum p counted 
from the K point for which a photoexcited e-h pair is sub- 
ject to relaxation by production of interband e-h pairs (carrier 
multiplication). The thick solid curve represents the contour 
Ep + e'lp — 2 eV, as determined by Eq. Q. The gray area 
shows the directions of p where both the electron with mo- 
mentum p and the hole with momentum — p are subject to 
relaxation; along the directions in the hatched area, only the 
hole is subject to relaxation. 



III. DISCUSSION 

Let us discuss some experimental implications of the 
obtained results. In an optical experiment, the incident 
photon of the frequency uj produces an electron with mo- 
mentum p and a hole with momentum —p. Their ener- 
gies satisfy + e'ip = uj, which constrains p to a trigo- 
nally warped circle, shown in Fig. [5] for cu = 2 eV. The 
direction of p is determined by the photon polarization. 
If the excitation density is low, one can neglect intra- 
band collisions between the photoexcited electrons and 
holes. According to the above results, when the band 
curvature is sufficiently large, the electron can relax by 
producing interband e-h pairs (carrier multiplication) if 
the direction of its momentum (the polar angle tp) satis- 
fies Cip = —Cs COS 3(y3 — Co > (shown in Fig. [s] by the gray 
area). The hole can relax if — Ctt-^ = — Cs cos3(/5-|-Co > 
(gray and hatched areas on Fig. [5]) . 

The discussed anisotropy of the electronic relax- 
ation has some implications for the two-phonon Raman 
scattering, who se intensity is suppressed by electronic 
relaxation.l^sEni^ Namely, it favors the electronic states 
near the KM direction (white sectors in Fig. [s] the states 
not subject to relaxation) to provide the dominant contri- 
bution to the two-phon on Ra man intensity, as has been 
observed experimentally.'^^J^ Another mechanism favor- 
ing the electronic states near the KM di rectio n is the 
anisotropy of the electron-phonon coupling.'^SlS! 

If the band curvature is too weak so that the quasi- 
particle state does not decay, its spectral weight Zp is 
still reduced, Zp < 1, since part of the spectral weight is 
transferred to three-particle excitations. It means that an 
initial excitation, produced by a short optical pulse, has 
a finite probability 1 — Zp to produce many-particle ex- 
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FIG. 6. (a) The basic elements of the diagrammatic tech- 
nique (the solid line represents the electronic Green's func- 
tion, the wavy line represents the Coulomb interaction) and 
the corresponding analytical expressions, (b) The two dia- 
grams (direct and exchange) contributing to ImEp(e) in the 
second order in e^/v. (c) Electronic polarization operator, 
(d) The RPA series for the self-energy. 

citations, the typical time of the processes being ~ l/7p 
(since 7p is the typical energy scale of the features in the 
spectral function). Thus, on average, the total number of 
e-h pairs per absorbed photon will exceed unity, so one 
can still speak about carrier multiplication even in the 
regime of weak band curvature. 
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Appendix A: Self-energy and polarization operator 
in the Dirac approximation and beyond 

1. General remarks about the calculation 

The calculation is performed using the standard zero- 
temperature diagrammatic technique whose basic ele- 
ments, the single-particle matrix Green's function [e — 
h{p) + iO+ signe]"^ with the 2x2 matrix h{p) given by 
Eq. ([3]), and the Coulomb interaction Vq = 27re^/|q|, are 
shown in Fig. |6]^a). 

The self-energy is also a 2 x 2 matrix. In the Dirac 
approximation, it can have components proportional to 
the unit matrix or to the scalar product p ■ cr, due to 
isotropy of the problem. Equivalently, the self-energy 
can be represented as 

2|p| 2|p| ^P^'^' 
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Pi 



FIG. 7. The initial electron momentum p (thick solid arrow) 
with its polar angle with respect to the KM direction (thin 
dotted line), and the final momenta of the two electrons pi , p2 
(thin solid arrows) and of the hole p — pi — P2 (thin dashed 
arrow) . 



where Sp(e),I]p(e) are scalar functions. The matrix co- 
efficients in front of them are the projectors on the two 
eigenstates of the Dirac Hamiltonian with momentum p. 
Beyond the Dirac approximation, the matrix structure 
of the self-energy becomes modified. However, this mod- 
ification represents a regular correction, proportional to 
^0, Csi so it is neglected in all calculations, since our pri- 
mary interest is the singularity at e — )■ vp. 

The singularities in the self-energy at e and in 

the polarization operator at a; — > come from nearly 
collinear processes, i. e., when all momenta are directed 
almost along the same line. Thus, all angular factors 
resulting from overlaps of eigenstates with different mo- 
menta, can be dropped, as they produce small regular 
corrections to the main singular behavior. This signifi- 
cantly simplifies the calculations. 



2. Second-order self-energy in the Dirac 
approximation 

Upon integration over the internal energy variables, 
the sum of the two diagrams in Fig. |6|^b) at e — >■ can 
be written as 



Imi;p(e) 



(27r)2 (27r)^ 



^-P2 I 



X 7r(5(e 



pi 



P2 



P-PI-P2; 



(Al) 



In the Dirac approximation = w|p|, so the (S-function 
constrains e/v to be equal to the sum of the lengths of 
the three thin arrows in Fig. [7] The triangle inequality 
ensures that this is possible only when e/v > p, the length 
of the long thick arrow. At e — >■ vp, the directions of 
Pij P2, P — Pi — P2 should approach the direction of p. 

Let xip,X2P be the projections of Pi,P2 on p, and 
yiP, y2P the projections on the orthogon al di rection. The 
main contribution to the integral in Eq. ( Al ) comes from 
the region xi,X2 > 0, a;i -I- a;2 < 1 and \yi\, \y2\ ^ 1. 



Then, we can approximate |p — Pi,2| ~ p(l — 2^1.2) in the 
Coulomb matrix elements, as they are nonsingular in the 
colhnear limit — > 0. In the energy J-function, yi_2 
are kept to the second order: 




which gives 



2 \ 2 

Iml]p(e)« - ( - ) ^ X 



1 1-2:1 
X J dxi J dx2 




N 



(l-a:i)2 (l-a;i)(l-a;2) 



X 



dyi dy2 S 



e-vp yl yl (?/i+y2)V2 



vp 



2xi 2X2 1 — Xi ~ X2 



The y-integration is performed using the general relation 

2n \ 2n 



°f / 2ri \ 2n 



(A2) 



\/detn(n- 1)!' 



valid for any positive-definite 2n x 2n matrix A. The 
remaining x-integration. 



ImEp(e) 



1— xi 



vp 



e — up) / dx\ \ dx2 X 



X ^a;ia;2(l ~ xi - X2) 



N 



(l-Xiy (1-Xi)(l-X2) 



(A3) 



straightforwardly gives Eqs. (6a), (6b I 



Eq. ( A2 ) also determines the suppression by energy- 



momentum conservation of higher-order contributions to 
Inil]p(e — > vp) corresponding to emission of n electron- 
hole pairs. Indeed, this is precisely the kind of integral 
one obtains for the perpendicular components of the 2n 
momenta in the collinear limit, with s oc {e — vp). 



3. Second-order self-energy beyond the Dirac 
approximation 

Neglecting the trigonal warping correction to the 
single-particle eigenstates, one can again use Eq. (Al) 
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with the quasiparticle dispersions from Eq. Q to cal- 
culate the self-energy. Using the same notation as in 
Appendix |A 2[ we expand the quasiparticle dispersions 
to the second order in yi,y2- Consider, for example, ep^: 

= vpxi - (Co + Ca cos 3ip)p'^xl + KzP^xiyi sin3(p + 

^ - CoP +2 CsP cos 3.^ I . 

(A4) 

The expression for has a similar structure. The 

terms linear in j/i, 7/2 can be removed by a shift of yi,y2- 
Neglecting the terms ~ (C3P^)^/('^p) ^ well as correc- 
tions to the determinant in Eq. (|A2|, we obtain the same 



Eq. ( A3 ) , but with a modified integration domain: in ad- 
dition to the conditions xi^X2 > 0, xi + X2 < 1, there is 
another condition 

e-vp- + 2^/(1 - xi)(l - X2) > 0, (A5) 

where = —Co ~ Ca cos Stp. Even though we could 
not evaluate the corresponding integral analytically, some 
general properties of it can be established. First, for suffi- 
ciently large negative e—vp < —\C^\p'^ the integration do- 
main is empty, so ImSp(e) = 0. Second, f or su fficiently 
large positive e~vp > |C<pb^i the condition (A5) becomes 
redundant, the integration domain coincides with that in 
Eq. (A3 1, so Iml]p(e) = 7p. Finally, the condition (A5) 



becomes its opposite upon simultaneous change of sign 
oi e — vp and Cip- Thus, if one defines a function Fn{z) 
as 



Dirac approximation in Sec. A 5 In the coUinear ap- 



proximation the angular factors can be neglected, so the 
polarization operator is given by 



nq(^) 



(27r)2 



N 



(A8) 



q-p 



Let us denote by xq the projection of p on q. and by yq 
on the orthogonal direction. Then, 



nq(^) 



Nq 



dx dyx 



y 

2\x 



y 



Nq 



2\l-x\ 
y^2\x{l - x)\dx 



vq 



1 



Nq-" 



2TTyj2vq{vq — w) 



\/ x{\ — x) dx ■ 



Nq r ^7/2 -1/Adu 

2nv J Ju\2u- uj/(vq)]' 
1/2 

The second term is nonsingular at cj — > vq, so it can be 
ignored (the divergence of the integral is spurious, being 
a consequence of the coUinear approximation). The first 
one gives -{Nq^/16)/^j2vq{vq ~ w), which is precisely 
Eq. (IA7| at w -> vq. 



Fn{z) = J In{xi,X2) dxidx2 

X J 6{z — 1 + 2(1 - xi){l — X2))lNixi,X2) dxidx2, 

(A6a) 

In{xi,X2) = d{xi) 6{x2) 6{l ~ xi- X2) X 

(A6b) 



N 



_(l-a;i)2 (l-a;i)(l-a;2) 
then ImEp(e) is given by Eq. ([s]). 



4. Polarization operator in the Dirac 
approximation 

The exact expression for the polarization operator in 
the Dirac approximation is known since long agd^ 



nq(w) 



iVg2/i6 



(A7) 



Still, here we give its simple derivation near the singu- 
larity at a; — > vq, which will be generalized beyond the 



Polarization operator beyond the Dirac 
approximation 



As in Sec. A 4 let us start from Eq. ( A8 1 and denote by 



xq the projection of p on q, and by yq on the orthogonal 
direction. As in Sec. |A 3[ let us expand the energies to the 
second order in y, see Eq. ( A4). In the interval < a; < 1, 



which determines the main singularity, we have 



q p 



vq+C,^q^{2x-l) + 



vq 



1 



1 



1 



(y-2/o)^ 



where y^ is an x-dependent shift. The integration over y 
is straightforward, the subsequent integral over x reduces 
to eUiptic integrals, to give Eqs. (|9a|), ([9b|). 



Let us mention some properties of the function V{z), 
defined in Eq. (9b). Many of them can be deduced di- 



rectly from the integral representation. 



V{z) 



du. 



(A9) 



For z < —1, 'P{z) is purely real, at 2 > —1 it acquires 
a positive imaginary part whose sign is fixed by the re- 
quirement of the analyticity in the upper half-plane of 
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the complex variable z, and V^z) is purely imaginary 
at z > 1. The real and imaginary parts are related by 
Re J'(-z) = Im J"(z). At z ^ -oo, V{z) 1/7^. At 
z — > ±1, V{z) has a weak singularity: 



V{z) 



8V2 V^(z + 1) 



In- 



A 



z + 1 



+ i7r6l(z + 1) 



(AlO) 

with A ^ 1. Finally, there are two integral relations, 
valid at ?; > 1: 

V 

Im'P(z) = TT, (Alia) 



/ 



dz ImP(z) 



Vy^z \V{zW 



Try 
2 ■ 



(Allb) 



The first relation can be proven by using Eq. ( A9 ) and 



interchanging the order of integration, while the second 
relation has been verified numerically. 



6. RPA self-energy in the Dirac approximation 

Upon integration over the internal frequency variable 
of the RPA diagrams, the self-energy becomes 



Imi;p(e) 



^ Im ■ 



(A12) 

The imaginary part of the dressed polarization operator 
is given by 



Im 



nq(^) 



1 - v^n^iio) 



16 ufi — xP-q^ 



- vq), 



= 1 - 



(A13a) 
(Al3b) 



\ 8 w 

Note that can become negative, which does not rep- 
resent any problem; the main effect of is to make the 
denominator nonzero at a; — > vq. Because of this, in the 
resulting expression for the self-energy, 



Imi;p(e) 



— \ (fqe{e - v\p - q| - vq) x 



V^(e-f|p-q|)2-i'2g^ 
(e — w|p — q|)^ — w^g^ ' 



(AM) 



one can simply set e = up in the denominator and re- 
place the latter by (u^ — v'^)q'^ provided that |e — vp\ <C 
{Ne^/vfvq. This condition will provide a lower cutoff 
for the integration over q. 

In the numerator, we approximate 



(e - wIp - q|)2 - v'^q^ « 2vq\\ (e - vp) - 



P - 1\\ 



where , q± are the components of q along p and per- 
pendicular to p, respectively, and \q±\ <C q\\. Then after 
the straightforward integration over q±, we obtain 



Im Sp(e) 



4(e — vp) f dq\\ jp - q\ 



ttN J (7|[ 
4:{e-vp) j^_P_ 



(A15) 



where the logarithmic divergence is cut off at 

f e — vp 

z;g„„, = max|e-z;p, ^^^y^^^ 

as discussed above. Thus, we arrive at Eq. Q. 

7. RPA self-energy beyond the Dirac 
approximation 



Let us again start from the general Eq. (A12). Beyond 



the Dirac approximation, we use Eq. (9a I and obtain 
T/q2nq(e - e^_q) 2^6^ VAlmV{z) 



Im 



1 



where we denoted 

.2N 2 



Ariz)\^ 

(A16a) 



A 



8v 

1 



vqii 



vp 



(A16b) 



vpql 



2911 (p- 911 

(Al6c) 

In the relevant region of energies, namely, where Eq. ([t]) 
is valid in the Dirac approximation, we have A^ I and 
z <C -4. Due to the latter condition and to the fact that 
Ini7'(z) ^ l/-\/z at z — > oo, one can neglect the unity in 
the denominator of Eq. ( Al6a), which then becomes 



Im 



V^n^ie-el_^) 
l-Fqnq(e-e: 



16 2\C^\v ImViz) 



p-q; 



TV 



(A17) 



Let us pass from the integration over q± to the one over z, 
which gives 



0{z„ 



x^/l — X dx X 



l)dz Im7'(z) 



y/zn 



\r{z)\^ 



vp - Cvpp^ii - ^y^ 



(A18a) 



(AlSb) 



The z-integration is performed using Eq. (Allb). On the 
mass shell, e = vp + Ci/^P^i the a;- integral in Eq. ( Al8a 
converges on the lower limit, and one arrives at Eq. ( |10[ ) 
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